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The quark behaviour in the background of intensive stochastic gluon field is studied. An approx- 
imate procedure for calculating the effective Hamiltonian is developed and the corresponding 

"ob" ' ground state within the Hartree-Fock-Bogolyubov approach is found. The comparative analysis 

of various model Hamiltonian is given and transition to the chiral limit in the Keldysh model is 

^ ! discused in detail. 

We study the quark (anti-quark) behaviour while being influenced by intensive stochastic gluon 
I field and work in the context of the Euclidean field theory. The corresponding Lagrangian density is 
^ the following 



£e = Q {iltiDy, + im) q , (1) 



here q (q) — are the quark (anti-quarks) fields with covariant derivative _D M = d^ — igA^t" 1 where A a ^ is 



the gluon field, t a = X a /2 are the generators of colour gauge group SU (N c ) and m is the current quark 
mass. As the model of stochastic gluon field we refer to the example of (anti-)instantons considering 
^h. an ensemble of these quasi-classical configurations. On our way to construct an effective theory 
(which usually encodes the predictions of a quantum field theory at low energies) the assumptions 
^ \ done are not of special importance. However, what is entirely restrictive to fix the effective action 
(Nj ; at really low energy (i.e. low cutoff) up to a few coupling constants is an idea to neglect all the 
^ ■ contributions coming from gluon fields A ex generated by the (anti-) quarks. 



A ex A . 



qq : Actually, it means the removal of corresponding cutoff(s) from consideration, but by the definition 
O of an effective theory this operation does not pose itself. Then the corresponding Hamiltonian 
>• ! description results from 

X'- H = nq-£ E , n = ^ = zq + , (2) 

dq 

■ ' and 

= -q («7V + im) q , (3) 

for noninteracting quarks. In Schrodinger representation the quark field evolution is determined by 
the equation for the quark probability amplitude $ as 

4? = -HV , (4) 

with the density of interaction Hamiltonian 

V s = q(x)t a ^A^t,x)q(x) . (5) 

The explicit dependence on "time" is present at the gluon field only. The creation and annihilation 
operators of quarks and anti-quarks a + ,a, b + ,b have no "time" dependence and consequently 

q a i(x) = J (2| p ]|) 1 /2 [ a ( p > s ' c ) Uai ( p i s ' c ) ^ + b+ (P' s > c ) Vai ( p i s > c ) e ~ iPX ] ■ ( 6 ) 
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The stochastic character of gluon field (which we supposed) allows us to develop the approximate 
description of the state \I/ if the following procedure of averaging 



/ dr *(r)/ t 
Jo 



is intoduced. With this procedure taken the futher step is to turn to the approach of constructing 

a density matrix (v]>\I/). However, here we believe that at calculating the ground state (or more 
generally with quasi-stationary state) it might be sufficiently informative to operate with the averaged 
amplitude directly. Then in the interaction representation \1/ = e H °'$ we have the equation for state 
$ as 

$ = -y$ , V = e Hot V s e- Hot . (7) 

Now the "time" dependence appears in quark operators as well and after averaging over the short- 
wavelength component one may obtain the following equation 

/•oo 

<$(*)> = + dr (V(t)V(t - r)> (*(*)> . (8) 



The limitations to have such a factorization validated are well known in the theory of stochastic 
differential equations (see, for example, pQ). The integration interval in Eq.(8) may be extended 
to the infinite "time" because of the rapid decrease (supposed) of the corresponding correlation 
function. Now we are allowed to deal with amplitude ($(£)) in the right hand side of Eq.(8) instead 
the amplitude with the shifted arguments in order to get an ordinary integro-differential equation. In 
the quantum field theory applications it is usually difficult to construct the correlation function in the 
most general form. However, if we are going to limit our interest by describing the long-wavelength 
quark component only then gluon field correlator (A a {x) A b (y)) may be factorized and as a result we 
have 

/* /*oo r 

($(*)) = Jdx q{x, t) t a la q(x, t) J o dr J dy q(y, t-r) t b lu q(y, t-r) g 2 (A«(t, x)A b u (t-r, y)) ($(t)> . 

Having assumed the correlation function rapidly decreasing in "time" we could ignore all the retarding 
effects in the quark operators. Turning back to the Schrodinger representation we have for the state 
amplitude x = e~ Hot (Q) the following equation 

X Hind X ) 

r roo , (9) 

n ind = -q (nV + im) q-q t a la q J dy q' t b lu q' j dr g 2 (A a a A b ) , 

with q = q(x), q = q(x), q' = q(y), q' = q(y) and A a a = A^(t,x), A' b = A b u (t — r,y). Now the 
correlation function might be presented as 



j o drg 2 (A a a A b ) = 5 ab F uu (x-y) 



with the corresponding formfactors F uu (x — y) = S uu I{x — y) + Jnu{ x ~y)- 111 our consideration we 
ignore the contribution of the second formfactor spanning on the components of the vector x — y. 
Thus, on output we receive the Hamiltonian of four-fermion interaction with the formfactor rooted in 
the presence of two quark currents in the points x and y. With this form of the effective Hamiltonian 
we could apply the Hartree-Fock-Bogolyubov method [2j to find its ground state as one constructed 
by the quark-anti-quark pairs with the oppositely directed momenta 

, 1 (10) 

a + (p, s, c) & + (-p, s, c) + a(p, s, c) b(-p, s, c) 
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where the parameter 9(p) characterizes the pairing strength. Introducing the creation and annihila- 
tion operators of quasi-particles A = T a T _1 , A + = T a + T _1 , B = T b T" 1 , B + = T b + T~ 1 , we 
can rewrite the quark (anti-quark) operators as 



q[x) 
q(x) 



dp 



(27r) 3 (2|p 4 |)V2 
dp 1 



(27T) 3 (2| P4 |) 1 /2 

with the quasi-particle spinors 



A{p, s, c) U{p, s, c) jP x + B+(p, s, c) V(p, s, c) e { P x 
A + (p, s, c) U(p, s, c) e-*P* + B(p, s, c) V(p, s, c) e^ x 



U{p,s,c) = cos y-J u(p, s, c) - sin ^-J v(-p,s,c) 
V(p,s,c) = sin J u(-p,s,c) + cos (~ J v(p,s,c) 



where U(p, s, c) = U + (p, s, c) 74, V(p, s, c) = V + (p, s, c) 74. Minimizing the mean energy functional 
one is able to determine the angle 9 magnitude 



d{a\H ind \a) 
d9 







(12) 



Dropping the calculation details out we present here the following result for the mean energy as a 
function of the 9 angle 

.2 

COS I 



I \TT I \ f d P 2N C P4 n 



-G 



(2tt)« 



(3/-J)t^ + (4/ -J)- Pg 



Tfl \ I Tfl \ 

sin 9 cos 9 sin 9' cos 9' ) + 

v J V 9 



\vAm\ 



n 771 ■ n\ nl 111 ■ nl 

cos y H sin # cos 9 H sin 

v P \ Q j 



(13) 



here the following designations are used p = |p|, g = |q|, / = I(p + <jr), Jy = Jij(p + q), J = Z)f=i ^ii> 
p 2 = q 2 = —m 2 , 9' = 9{q) where G is the constant of corresponding four-fermion interaction (the 
relevant details can be found in [3J). The first integral in Eq. ( 1T31) comes from free Hamiltonian, and 
we make a natural subtraction (adding the unit) in order to have zero mean free energy when the 
angle of pairing is trivial. 



Nambu— Jona-Lasinio model 

In order to get an idea of the parameter scales we continue with handling the model in which 
the formfactor behaves in the coordinate space as I(x — y) — S(x — y), J^ u = 0, dropping contribu- 
tion spanned on the piqj tensor also. Actually, it corresponds to the Nambu- Jona-Lasinio model [I]. 
As well known the model with such a formfactor requires the regularization and, hence, the cutoff 
parameter A comes to the play 



W 



dp 



\P4 



COS! 



G-r^-r (SHI (9 



m 



n\ f A d <l Q ( ■ nl 
cos 9 I - — -, — r sin 9 

p J J (2vr) 3 \q 4 \ \ 



m 



cos 9' 



■ (14) 



We adjust the NJL model with the parameter set given by Hatsuda and Kunihiro [4 J in which 
A = 631MeV, m = 5.5MeV. One curious point of this model is that the solution for optimal angle 9 
in the whole interval p G [0, A] can be found by solving the simple trigonometrical equation 



(p 2 + m 2 ) sin 9 — M q (p cos 9 + m sin 9) = 



(15) 
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Figure 1: Phase portrait of the Keldysh model, sin 9 as a function of momentum p(MeV). The dotted 
curve corresponds to the solution with the negative values of angle in the chiral limit m = 0. 



with the dynamical quark mass 



M„ 



f A dp p ( m 

2G / t t"tt i r SlIlS COS 9 

(2tt) 3 |p 4 | ^ p 



(16) 



Eventually the results obtained look like M q = —335 MeV for dynamical quark mass and (&\qq\<j) 
—i (245 MeV) 3 for the quark condensate with the following definition of the quark condensate 



{cr\qq\(T) = — r 

71 



/ dp - — | (p sin 9 — m cos 9) . 
Jo \p 4 \ 



(17) 



The Keldysh model 



Now we are going to analyse the limit, in some extent, opposite to the NJL model, i.e. we are 
dealing with the formfactor behaving as a delta function but in the momentum space (analogously 
the Keldysh model, well known in the physics of condensed matter [5]), I(p) = (2tt) 3 S(p). Here the 
mean energy functional has the following form 



W(m) 



dp 



\p 4 \ (l-cos0)-G 



P 



\Pa\ 2 



sin 9 cos 9 

k P j 



contrary to the NJL model there is no need to introduce any cut off. The equation for calculating 
the optimal angle 9 becomes the transcendental one 



\p 4 1 3 sin 9 — 2G (p cos 9 + m sin 9) (p sin 9 — m cos 9) = , 



(19) 



and, clearly, it is rather difficult to get its solution in a general form. Fortunately, it is much easier 
and quite informative to analyse the model in the chiral limit m = 0. There exist one trivial solution 
9 = and two nontrivial ones (for the positive and negative angles) which obey the equation 



a p 
COs9= 2G 



(20) 



Obviously, these solutions are reasonable if the momentum is limited by p < 2G. Then for the mean 
energy we have ^±(0) = if the quark condensate defined as (a\qq\a)(0) = — rf^ — . For 
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the trivial solution the mean energy equals to zero together with the quark condensate W'o(O) = 0, 

(a\qq\a) o(0) = 0. Introducing the practical designation sin^ = — "^2x1/2 which characterizes 

(p + M e ) 

the pairing strength by the parameter Mg we have, for example, for the nontrivial solution Mg = 
(4G 2 — p 2 ) 1 ^ 2 - In order to compare the results with the NJL model we fixed the value of four-fermion 
interaction constant as M e (0) = 2G = 335 MeV. It is interesting to notice that the respective energy 

becomes constant E(p) = Jp 2 + M|, E{p) = 2G. 



t — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 j_ 1 — 1 — 1 — 1 — 1 — 1 — 1 — r 




Ptf P 



Figure 2: The optimal angle 6 as a function of momentum p(MeV). The solid line corresponds to 
the NJL model and the dashed one to the Keldysh model. The current quark mass is m = 5.5 MeV 
and p e ~ 40 MeV. 

After having done the analysis in the chiral limit which is shown by the dotted line in Fig.l we 
would like to comment the situation beyond this limit. The evolution of corresponding branches 
is available on the same plot [TJ One solution denoted by A is developing in the local vicinity of 
coordinate origin and for small values of quark mass this domain is practically indistinguishable. 
In order to make it noticeable (to have a reasonable resolution on the plot) the quark mass was 
put as m = 50 MeV. Besides, there are two solutions a and b in the domain denoted by /, three 
solutions denoted by 1, 2, 3 in the domain II and one solution B in the domain III. The minimum 
of mean energy functional can be realized with the piecewise continuous functions. At the local 
vicinity of coordinate origin we start with the solution branch A, then relevant solution passes to 
the branch a or b interchanging its position from a to b in any subinterval. But in any case there 
is only one way to continue the solution at streaming to the infinite limit and it is related with the 
branch B where the angle is going to the zero value. As to the functional ( TTBl) the contribution 
of the term proportional to the cosine in the second parenthesis is divergent even if the angle 9 is 
zero. It means the mean energy out of chiral limit goes to an infinity at any nonzero value of quark 
mass. The same conclusion is valid for the chiral condensate. In principle this functional could be 
regularized and corresponding continuation might be done but it is out of this presentation scope. It 
is not difficult to demonstrate the similar discontinuities of functional are present, for example, for 
Gaussian I{x) = G exp (—a 2 x 2 ), and exponential I{x) = G exp (—a \x\), formfactors and they are 
present even in the NJL model but this fact is masked by the cut off parameter. 

Comparing the optimal angles in the NJL and Keldysh models (see Fig. 2) it is interesting 
to notice that the formation of quasiparticles becomes significant at some momentum value close 
to the origin pg ~ 40 MeV but not directly at the zero value. It is clear the inverse value of 
this parameter determines the characteristic size of quasiparticle. Parameter Mg as a function of 
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Figure 3: The parameter M 9 (MeV) as a function of momentum p(GeV) corresponding to the best 
fit of the NJL data M q = 335 MeV, (cr\qq\a) = -i (245 MeV) 3 . The solid line corresponds to the 
Gaussian formfactor in chiral limit and the dashed line corresponds to the magnitude of current 
quark mass m = 5.5 MeV. The exponentially behaving formfactor is represented by the dotted lines 
and p ~ 150 MeV. 

momentum p corresponding to the best fit to the NJL data M q = 335 MeV, (a\qq\a) = —i (245 
MeV) 3 is shown in Fig. 3. The solid line corresponds to the Gaussian formfactor in the chiral limit 
and the dashed one shows the same dependence for the current quark mass m = 5.5 MeV. This 
dependence for exponential behaviour of formfactor is presented by the dotted lines on the same plot 
(the characteristic angle is pg ~ 150 MeV in this case). Analysing the discontinuity of mean energy 
functional and quark condensate we face some troubles at fitting the quark condensate, for example. 
However, the dynamical quark mass and quark condensate are nonobservable quantities and it is 
curious to remark here that although the mean energy of the quark system is minus infinity the 
meson observables are finite and even in Keldysh model the mesons are recognizable with reasonable 
scale and we can in principle make a fit for this observables. 
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